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ABSTRACT 

We present a formalism to infer the presence of merging activity by comparing the observed time derivative 
of the galaxy stellar mass function (MF) to the change of the MF expected from the star formation rate (SFR) 
in galaxies as a function of galaxy mass and time. We present the SFR in the Fors Deep Field as a function 
of stellar mass and time spanning 9 < logM* < 12 and < z < 5. We show that at z > 3 the average SFR, 
M*, is a power law of stellar mass (M* oc M° 6 ). The average SFR in the most massive objects at this redshift 
are 100 - 500 M Q yr _1 . At z ~ 3, the SFR starts to drop at the high mass end. As redshift decreases further, 
the SFR drop at progressively lower masses (downsizing), dropping most rapidly for high mass (logM* > 11) 
galaxies. The mass at which the SFR starts to deviate from the power-law form (break mass) progresses 
smoothly from \ogM\ > 13 at z ~ 5 to \ogM\ ~ 10.9 at z ~ 0.5. The break mass is evolving with redshift 
according to M\{z) = 2.7 x 10 10 (1 +z) 2A . We directly observe a relationship between star formation history 
(SFH) and galaxy mass. More massive galaxies have steeper and earlier onsets of star formation, their peak 
SFR occurs earlier and is higher, and the following exponential decay has a shorter e-folding time. The SFR 
observed in high mass galaxies at z ~ 4 is sufficient to explain their rapid increase in number density. Within 
large uncertainties, we find that at most 0.8 effective major mergers per Gyr are consistent with the data at 
high z, yet enough to transform most high mass objects into ellipticals contemporaneously with their major star 
formation episode. In contrast, at z < 1.5 and at M* > 11, mergers contribute 0.1 - 0.2 Gyr -1 to the relative 
increase in number density. This corresponds to ~ 1 major merger per massive object at 1.5 > z > 0. At 
10 < logM* < 11, we find that these galaxies are being preferably destroyed in mergers at early times, while 
at later times the change in their numbers turns positive. This is an indication of the top-down buildup of the 
red sequence suggested by other recent observations. 

Subject headings: surveys — cosmology: observations — galaxies: mass function — galaxies: evolution — 
galaxies: star formation 



1. INTRODUCTION 

Galaxy formation and evolution in a cold dark matter 
(CDM) dominated universe can be described in the most gen- 
eral terms as being driven by the gravitationally induced hi- 
erarchical growth of structure in the dark matter component 
coupled to the baryonic matter within the merging and ac- 
creting dark matter halos. While the dark matter is collision- 
less, baryonic matter dissipates energy by radiative cooling 
and thereby is able to form stars. Feedback effects are thought 
to self -regulate the latter process. 

Directly observing the fundamentally important hierarchi- 
cal growth of dark matter halos is not possible at the present 
time. Information about the growth of the halos has to be de- 
duced from observations of galaxies or clusters. Large gravi- 
tational lensing surveys detecting cosmic shear will allow the 
study of dark matter structure formation directly in the future. 
In the meantime, it has proven difficult to deduce halo growth 
from observations of galaxy growth. 

Galaxies themselves grow in stellar mass firstly by star for- 
mation, and secondly by accreting smaller galaxies and by 
merging with similarly sized galaxies. Both of these processes 
need to be mapped out if we wish to fully understand the as- 
sembly history of galaxies and their host dark matter halos. 

Progress has been made in recent years in measuring the 



build-up of ste llar mass density from redshift z ~ 6 to the 
present epoch fBrinchmann & Elli s] l2000t iDrorv et all 1 20011 
I Cohenl 12001 [Dickinson et all 120031: iFontana et al J booa 
IRudnick et al. 2003; Glazebrook et al. 2004; Chapma n et all 
I 2005J IRudnick et ail l2006t lEyles et all 120071: iGrazian et all 
120071 IStark et all 120071k Also, the stellar mass function 
has been mapped locally (ICole et all 1200 It [Bell et alJl200l 
iPerez-Gonzalez et alJ l2003a|) and in some detail to z ~ 1 .5 



(IDrorv et all l2004t iFontana et alJ 1 20041 ; [Bundv et al. 
Borc het all 120061; Pannella et alJ \200a Bundv et al 



2005 



2006; 



iPozzetti et alJl2007l) with generally good agreement between 
different datasets. To some lesser detail and ac curacy deep 



surveys have provided data spanning < z < 5 dDrorv et al 
2005b IConselice et alJ 120051 IFontana et al. 20061: lYan etal 



20061: IGrazian et alJl2007l) and even some estimates at z 
dBouwens et al.ll2006l) . These mass functions tend to agree 
fairly well at z < 2, while the integrated mass densities may 
differ by a factor of ~ 2 at higher z (see, e .g., comparisons in 
IDrorv et al. 2001 and IFon tana et al. 2006]). In summary, data 
are now available that describe the mass in stars from an early 
epoch where only ~ 1% of the current mass density was in 
place until the present time. 

Deep optical and near-infrared surveys, as well as 
mid-infrared data from the Spitzer Space Telescope and 
sub-mm surveys have greatly expanded our understand- 
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ing of the cosmic star forma tion (SF) hi s tory. The in- 
tegrated star formation rate dMadau et al.l 1 19961) is now 
quite well constrained to z ~ 4 (with some estimates at 
even higher z) from a multitude of studies employing dif- 
ferent methods at wavelengths from the UV to the far- 



IR dLillv et al. 


19961 iMadau et all Il998t iFlores et alJ |1999t 


Pei et al. 1999 


Kauffmann et al.ll2003al iHeavens et al.ll20o4 


Gabasch et al. 2004b; Chapman et al. 2005; Bell et alj|2005l; 
Wang et al. 2006; Thompson et al. 2006; Hopkins & Beacom 


2006 
2006 


, iNagamine et al.ll2006l lYan et al.ll2006l iBouwens et all 



In the recent literature, more attention is being de- 
voted to the star formation rate (SFR) as a function 
of stellar mass, usually in the form of the specific star 
formation rate (SSFR) as a function of stellar mass 
dCowie et al.lll996t iBrinchmann & Ellisl 120001: iFontana et al 



2003b; 



Jrinchmann et al 
z.u05t ISalimetal 



2003; Perez-Gonzalez et al. 
2004[ lBaueretal.1 1 20051: lJuneauetal 
20051 (Feulner et al l l2005allbl: iPerez-Gonzalez et alJ 120051; 
Feulner et al.ll2006tlBell et1dl2007l:IZneng et al. 2007|). 

Progress on the merger rate and its evolution with 
redshift has been somewhat slower. Observational re- 
sults obtained from galaxy pair statistics vary consid- 
erably, finding merger rates oc ( 1 + z ) m with < m 
4 in the redshift range < z < 1 dZepf&Kool 119891; 



Burkev et al. 1994; Yee & Ellingson 1995; Neuschaefer et al 
19971; iLe Fevre et alJl2000t ICarlberg et a l. 2000; Patton et al 
2002[ [Lin et al l 120041; iBundv et all l2004t iDe Propris et a~ 
2005] IConselicel 120061: iBell et al.l 12006M: |Kartaltepe et al 
20071) . Early studies tend to give somewhat higher numbers 
than the more recent work. The diversity of these obser- 
vational merger rates is likely to be due to differing obser- 
vational techniques, differing criteria for defining a merging 
pair, and survey selection effects. To this date, no results have 
been published for significantly higher redshifts. Note that 
Cold Dark Matter (CDM) N-body simulati ons predict halo 
merger rates oc (1 +z) m with 2.5 < m < 3.5 (iGovernato et al.1 
1999: Gottlober et al]l2TMl2002l) . 

Little is known about the dependence of the merger r ate on 
galaxy luminosity or mass. Recently, IXu et al.l (120041) find 
evidence from the A" s -band luminosity function of galaxies in 
major-merger pairs that the local merger rate decreases with 
galaxy mass. 

In this paper, we present a method to infer the presence 
of merging activity and as such the underlying hierarchical 
growth of structure from a comparison of observations of 
galaxy mass functions with observations of star formation 
rates in the galaxies. Essentially, we compare a measure- 
ment of the time derivative of the mass function (the num- 
ber density of galaxies per unit stellar mass; (/>(M*,?)) with an 
expectation derived from time-integrating the (average) star 
formation rate of galaxies as a function of mass, M*(M*,t). 
The difference between the two can then be attributed to mass 
function evolution due to "external causes", which we iden- 
tify with accretion and merging (after taking the survey selec- 
tion effects into account). We present our results with caution, 
given the limitations of the current dataset. 

This paper is organized as follows: In § [2] we establish a 
formalism that compares the change in the stellar mass func- 
tion expected from the observed star formation rate as a func- 
tion of mass and time to the observed time derivative of the 
mass function. The aim is to deduce a merger rate (or, more 
precisely, the change in number density due to merging) as a 



function of galaxy mass and time from the difference of the 
two. We will call this quantity the assembly rate. § [3] and 
§ |4]present the mass function data and the star formation rate 
data, respectively. Also, suitable parametrizations of the data 
and its evolution with time are presented and discussed. § [5] 
discusses the total time derivative of the mass function. § [6] 
discusses the change of the mass function that is expected due 
to the star formation rate only. § [7] finally combines the for- 
mer two results to discuss the change in number density due 
to merging we infer from our data. All these are summarized 
and discussed in context in §[8] 

Throughout this work we assume VIm = 0.3, fl\ = 0.7, Hq = 
70kms _1 Mpc" 1 . 

2. DESCRIBING THE EVOLUTION OF GALAXY MASS 

A general description of the evolution of galaxy stellar mass 
with time can be given in terms of (1) internal evolution due 
to star formation in the galaxy, and (2) externally-driven evo- 
lution due to the infall of other smaller galaxies (accretion) 
or the merging with another galaxy of comparable size. The 
former can be measured by finding the average star formation 
rate as a function of stellar mass and time, M*(M* , t). The lat- 
ter cannot be measured reliably at this point. However, if one 
knew the time derivative of the galaxy stellar mass function 
(GSMF), dcj){M*,t)/dt, then it could in principle be derived 
by comparing the internal evolution due to star formation to 
the total observed evolution, and thereby inferred from the 
difference between the two. 

First, we note that the evolution of the GSMF due to star 
formation in the galaxies by definition preserves the number 
of galaxies. Hence the change of the GSMF due to star for- 
mation is described by a continuity equation in <j)(M* , t): 



d(j)(M,,t) 



dt 



SF 



d 

"dM, 



[</>(M*,t)M*(M*,t)] 



(1) 



Here, somewhat unconventionally, we have taken |sf to 
mean "due to star formation". 

This equation becomes intuitively clear if one considers that 
the continuity equation above means that the change in num- 
bers in a given stellar mass bin [M,M + dM) is given by the 
difference across this bin of the number of objects entering 
the bin from below and the number of galaxies leaving the 
bin to higher masses, both owing to their star formation rate. 
Hence the derivative with respect to mass of the number den- 
sity times the star formation rate (which gives the difference 
of mass flux across the bin) on the right hand side. 

The change with time of the GSMF, d(j)(M*,t)/dt, and 
the average SFR as a function of galaxy mass and time, 
M^M*,?), are in principle directly observable. We can there- 
fore infer the change in the GSMF due to all processes other 
than star formation (which we will identify with merging) 
from these two functions. We will call the resulting quantity 
the assembly rate: 



d(p(M*,t) 



dt 



d(KM*,t) d<p(M*,t) 



dt 



Hence, 



dt 



1 d 



Merge 



dt 



><9M* 



dt 



SF 



(2) 



(3) 



where we have left out the arguments for brevity and di- 
vided by <fi to obtain an equation for the relative change. 
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FIG. 1. — The stellar mass function (MF) as a function of redshift. Data are from the Fors Deep Field iDrory et al. 2005). Schechter-fits to the data are over 
plotted as solid lines. Confidence regions are indicated by shading (random uncertainties only). The low-z MF is repeated in the higher-z panels for comparison. 



Note that (tidal) stripping of stars from a galaxy (e.g. in 
clusters or in interactions) contributes to eq. [3] as well. As- 
suming that this loss of stellar mass is small compared to the 
overall growth rate of a galaxy, it can be neglected. 

The assembly rate is given in terms of the change in number 
density per unit time rather than in number of discrete merg- 
ing events per object per unit time (and per unit merger mass 
ratio). Therefore, we can only speak of an effective merger 
rate, equivalent in its effect on the number density to an equal 
number of same-mass mergers per unit time interval. An al- 
ternative is to think about this quantity as simply representing 
the average effective mass accretion rate. We cannot discern 
the mass spectrum of accretion and merging events that ac- 
tually cause the observed change in the mass function. It is 
also important to point out that a value of zero does not imply 
that no accretion and merging are happening, but only that the 
mass flux entering a mass bin through merging and the mass 
flux leaving that same mass bin through merging are equal, 
and hence the number density at this mass does not change. 
Values below zero imply that the number density decreases 
due to merging activity, and hence objects are preferably be- 
ing accreted onto more massive objects and thus destroyed. 
To emphasize the difference between what we measure and 
the merger rate expressed as the number of merger events per 
unit time and object, we will avoid using the term merger rate 
and prefer the name assembly rate. 

In the following sections we shall compile the necessary 
terms from GSMF and star formation data available in the 
literature. We will concentrate on the methodology and defer 
discussing selection effects and other deficiencies in the data 
to a later section where we proceed to plot eq. (0 and present 
our overall results. 

3. THE GALAXY STELLAR MASS FUNCTION 

The time derivative of the galaxy stellar mass function can 
be obtained in a fairly straight-forward way from observa- 
tional studies of the GSMF. Here, we use the stellar mass 



function data spanning < z < 5 by lDroryetaD (120051) from 
the Fors Deep Field (FDF). 

The FDF photometric catalog is published in iHeidt et all 
(2003) and we use the I-band selected subsample cover- 
ing the deepest cent ral region of the field as described in 
Gabas ciTet al.l d2004al) . This catalog lists 5557 galaxies down 
to / ~ 26.8. The latter paper also shows that this catalog 
misses at most 10% of the objects found in u ltra-deep ri- 
band observations found by Labb e et ail (120031) . Distances 
are estimated by photometr ic redshifts ca librated against 
spectroscopy up to z ~ 5 dNoll et al.l 120041) providing an 
accuracy of Az/(z, t/ „, ,.+ 1) < 0.03 with only ^ 1% outliers 
dGabasch et alj2004al) . The spectroscopy samples objects sta- 
tistically down to / = 25. 

The stellar masses are computed by fitting models of com- 
posite stellar populations of varying star formation history, 
age, metallicity, burst fractions, and dust content to UB- 
gRIZJK multicolor photometry. We assume a Salpeter stellar 
initial mass function. At faint magnitudes, the uncertainties 
in the photometric distances are dominated by the large un- 
certainties in the colors. These are taken into account in the 
fitting procedure and hence reflected in the uncertainties in 
the individual masses entering the mass function. We refer 
the reader to the original publications for full details. 

To lessen the impact of noise in the data on the time deriva- 
tives we are interested in, we choose to parameterize the 
evolution of the mass function, which we assume to be of 
Schechter form, by the following relations: 



<*M)m = „ [jL 



exp 



(4) 

M* I M* 



and 



0*(z) = pdl + z) p \ 

logM*(z) = qi +q 2 M\+z), 

a(z) = ao . 



(5) 
(6) 
(7) 
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FIG. 2. — Evolution with redshift of the parameters </>* and M* of the 
Schechter functions fit to the data shown in Fig.[T] The shaded areas show the 
1-cr confidence region of the fit as determined from Monte Carlo simulations 
(see text). Note that the derivatives w.r.t. z have much smaller uncertainties 
than the values of <j>*(z) and M*(z) themselves. 

First, we fit in dividual Schechter functions to the data from 
iDrory et all (I2005I) in each redshift bin (Fig. B). Then we fit 
eqs. (5)-(7) to the Schechter parameters as a function of red- 
shift, 4>*{z) and logM*(z), to obtain a smooth description of 
the redshift evolution of the stellar mass function (Fig. 
The error bars of each data point are used to weight the data 
during fitting. We determine confidence intervals for <fi*(z) 
and logM*(z) by Monte Carlo simulations. These confidence 
regions are indicated by shaded regions in Fig. [2] A conse- 
quence of parameterizing 4>*(z) and logM*(z) in this way is 
that in spite of the large uncertainties of these quantities, the 
constraints on their derivatives with respect to time are a lot 
tighter. 

In the redshift range < z < 2, where the mass function 
data are deep enough to constrain the faint-end slope of the 
mass function, a, it turns out to be consistent with being con- 
stant at a = —1.3. At higher redshift, the data do not allow to 
constrain the faint-end slope any m ore. This is consistent with 
the results by Fontana et al. 2 0061 

We therefore fix alpha at its low-z value for all redshifts. 
The consequences of this will be discussed below. 

The mass function and the Schechter-fits to the data are 
shown in Fig.[TJ and the best-fitting Schechter parameters are 
listed in Table [TJ The fit to the evolution of the mass func- 
tion with redshift as parameterized by eqs. (5)-(7) is shown 



TABLE 1 

Schechter fit parameters to the mass 
function 



z 


<t>* 

feLMpc 3 dex 1 


logM* 
ft 70 M © 


a 


0.50 


(2.01 ±0.10) x 10- 3 


11.25 ±0.02 


-1.3 


1.00 


(1.60 ±0.12) x 10~ 3 


11.22 ±0.03 


-1.3 


1.50 


(1.45 ±0.14) x 10~ 3 


11.16±0.03 


-1.3 


2.00 


(9.85 ± 1.12) x 10- 4 


11. 00 ±0.04 


-1.3 


2.75 


(9.70 ± 1.32) x 10- 4 


11.05 ±0.04 


-1.3 


3.50 


(6.50 ± 1.50) x 10-* 


10.99 ±0.05 


-1.3 


4.50 


(2.61 ± 1.71) x 10~ 4 


10.97 ±0.07 


-1.3 



in Fig. 12 The best fitting parametrization is 

4>*(z) = (0.0031 ±0.0002)(l+z)- 107±013 , 
logM*(z) = ( 1 1. 35 ±0.04)- (0.22 ± 0.03) ln(l+z), 
a(z) = -1.3, 

and we will use this form in the following. 

Using this best-fitting parametrization of the evolution of 
the mass function, we can extrapolate the values of the pa- 
rameters to z = 0. We obtain 0*(O) = 0.0031 /^ n M pc~ 3 dex~' 
which is consisten t with the value s obtained by ICole et alJ 
(12001b (0.003) and lBell et all d2003l) (0.0034) when these are 
rescaled to the same initial mass function and cosmology. The 
value we obtain for the characteristic mass, logM*(0) = 1 1.35, 
is higher than the values obtained in the mentioned literature 
(11.16 and 11.01, respectively); however, our value of the 
faint-e nd slope, a, is diff erent (-1.3 i n this work as opp osed to 
-1.1 in lCole et al.ll200Tl and -1.18 in lBell et al.ll2003h . If one 
takes the coupling between M* and a into account, a smaller 
value of a is degenerate with a higher value of M*, and there- 
fore our values agree very well with the literature at zero red- 
shift, once "corrected" to the same a along the lines of the 
degeneracy. 

In any case, the focus of this paper is to establish a method 
of disentangling star formation and merging contributions to 
galaxy stellar mass assembly. We shall therefore not be con- 
cerned too much about peculiarities of the mass function, in- 
stea d we refer the interested reader to the relevant literature 
(e.g.lBell et alJl2003l:lDrorv et al.ll2005l:IConselice et al.ll2005b 
Fonta na et alj|2006l and references therein). 

4. THE AVERAGE STAR FORMATION RATE 

To find the (average) star formation rate at each stellar 
mass and time, we t ake advantage of the data published by 
iFeulner et alJ (12005a) in the Fors Deep Field. 

We use dust-corrected UV-continuum emission around 
2800A to estimate star formation rates. The reason for doing 
so is twofold. First, there are no other homogeneous datasets 
that provide more robust star formation rates over the whole 
redshift range of interest. Second, the UV continuum lumi- 
nosity can be measured for every galaxy that has a stellar mass 
assigned to it. Therefore we have star formation rates for the 
exact same sample that is used to measure the stellar mass 
function. This at least provides some assurance of internal 
consistency. 

The dust corrections are estimated from the stellar popula- 
tion model fi ts to the multicolor p hotometry (see details and 
discussion in Feulner et al. 20 05ah . 

How reli a ble a re the UV-based star formation rates? 
iDaddi et ail d2007l) analyze multi-wavelength data (optical- 
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FIG. 3. — The upper plot shows the average star formation rate as a function of stellar mass in seven redshift bins ranging 0.5 < z < 4.5. Each panel shows the 
data of one redshift bin (filled circles). The confidence regions are indicated by shading. Fitting functions to the data (power laws with exponential cutoff at high 
mass; see text) are plotted as dashed lines. Each panel also shows the fitting functions in all other redshift bins as dotted lines for comparison. The lower plot 
shows the same data plotted against redshift. Each panel shows the star formation rate for a different stellar mass as a thick solid line. The dotted lines indicate 
the data for the other masses for comparison. Confidence regions are shaded. 
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TABLE 2 

Fit parameters to the star formation rate 

(eq.hj 



1 
1 

j1 

- 1 
± 1 



2.0 
1.8 
1.6 
1.4 

1.2E- 
1.0 

0.8 h . 

0.6 



" " I I " 

M* 1 = 2.66E+10(1+z) 2 - 1 



z 




log Mi 


P 




M yr 1 


M Q 




0.50 


13.1 ±1.6 


10.77 ±0.08 


0.69 ±0.11 


1.00 


20.1 ±1.0 


11.04 ±0.04 


0.56 ±0.08 


1.50 


43.6 ±2.8 


11.35 ±0.09 


0.52 ±0.07 


2.00 


88.2±3.1 


11.41 ±0.11 


0.61 ±0.07 


2.75 


132 ±3.9 


11.67±0.13 


0.59 ±0.05 


3.50 


336±12.7 


11.78(-0.25 + oo)' 


0.63 ±0.02 


4.50 


576 ±26.7 


12.95(-0.31 + oo)' 


0.66 ±0.01 



Values are unconstrained on the high side. They are 
consistent with a single power law with no break at all. 

NIR, Spitzer mid-IR, an d radio) at 1.4 < z < 2.5 in the 
GOODS region (see also iDaddi etaT]|2005l) . They conclude 
that the SFR can be estimated consistently from the UV light 
using SED or color-based dust corrections, with the UV-based 
estimates being within a factor of ~ 2 of the 24 /im-based es- 
timate for most galaxies (Sub-mm galaxies are an exception, 
as for those, the UV underestimates the SFR by at least an 
order of magnitude, but these are likely missed by the I-band 
selected sample used here). At high IR luminosities the dif- 
ference is higher, but it is often unclear whether it is due to 
an IR excess due to an AGN, which is likely because the UV 
estimate is in agreement with the radio-based estimate of the 
SFR. T hese conclus i ons ar e also supported by a similar anal- 
ysis bv lReddv et al.l d2005l) . These data only probe relatively 
massive galaxies, though. 

To estimate errors (or rather the range of possible values) 
of the average star formation rate as a function of mass, we 
proceed as following: instead of applying a further upward 
correction to all our star formation rates, we set the likely up- 
per envelope of the average star formation rate as the aver- 
age value plus the width of the distribution above the average 
(using the 68-percentile). The reasoning behind this is that 
the brightest UV sources at each mass are likely the least ex- 
tinct ones, and therefore they mark the upper limit. The likely 
lower confidence limit on the average star formation rate is 
calculated by estimating the maximum amount by which the 
average might be lowered due to low-SFR galaxies missing in 
the sample. To do this we need to assume a distribution of star 
formation rates in objects below our detection limit. We as- 
sume a star formation rate function power- law s lope o f - 1 .45 
following the only measurement by iBell et al.l d2007l) . Note, 
though, that this measurement is at z < 1 . We normalize to 
our mean densities and assume that the undetected population 
follows a power-law distribution of star formation rate of that 
slope at any given constant stellar mass and redshift. We then 
integrate the star formation rate function up to the point at 
which the objects would reach the detection limit of the FDF. 
The amount by which this missed star formation lowers the 
average value at each mass and redshift then defines the lower 
confidence limit. The region of confidence is shaded in Fig. [3] 

To reduce the influence of noise in the data on the deriva- 
tives we are interested in, we will parameterize the average 
star formation rate as a function of mass by an analytic ex- 
pression, just as we have done with the mass function. We 
choose a power law of the stellar mass with an exponential 
cutoff at high mass: 



FIG. 4. — Evolution of the characteristic mass, M' , above which the star 
formation rate as a function of mass begins to drop. 

The data along with this fit are shown in Fig. [3] The best- 
fitting parameters are listed in Table [2] Since such data on 
the star formation rate as a function of stellar mass and time 
have not been shown before in this form in the literature, we 
will take the time to digress from our goal of determining the 
assembly rate and discuss the star formation rate as it provides 
valuable insights. 

A few regularities in the data are apparent. At z > 3 the 
star formation rate monotonically rises with stellar mass to 
the highest masses we probe (logM, ~ 11.5). It is consis- 
tent with a single power law with exponent (3 ~ 0.6. The star 
formation rates in the most massive objects at this high red- 
shift are of the order of 100 - 500 Moyr" 1 , suggesting that 
these are the progenitors of modern day massive ellipticals 
seen during their major epoch of star formation. The average 
star formation rate drops most rapidly with redshift for high 
mass (logM* > 11) galaxies. Their period of active growth 
by star formation seems to be over by z ~ 1.5, and the star 
formation rate as a function of stellar mass turns over from in- 
creasing with mass to rapidly decreasing with both mass and 
time in objects with logM* > 11. 

At redshift z ~ 3, the star formation rates begin to drop at 
the high mass end while still being an increasing function 
of mass at lower masses. As redshift decreases further, the 
star formation rates drop at progressi vely lower masses. This 
behavior, usually termed downsizing (ICowie et al. 19%|), has 
been observed in numerous recent high-redshift studies (e.g . 
[ Bauer et al.ll2005l: Ivan der Wei et al.H"2005l:lBundv et al.ll2006l: 
Noes ke et al.1 l2007ab and is also seen in studies of the fos- 
sil r ecord of the stellar populations in the local universe 
(e.g. iHeavens et al.ll2004t (Thomas et al.ll2005l) . The mass at 
which the star formation rate starts to deviate from the power 
law form progresses smoothly from logM' ~ 13 at z ~ 5 to 
logMj ~ 10.5 at z ~ . 5. For comparison, in the local uni- 
verse IKauffman n et al.1 (l2003bl) find a characteristic mass of 
logM, ~ 10.4 at which the galaxy population changes from 
actively star forming at lower stellar mass to typically quies- 
cent at higher masses. This transition mass shows a smooth 
evolution to higher masses at earlier times up to very high red- 
shift (at > 3 the data are consistent with a single power law 
and no break mass at all). 

How does this "break mass" - the mass above which the 
average star formation rate begins to decline - evolve with 
redshift? Fig.|4]shows the evolution of the parameter M\ from 
eq. dHJ against redshift, as well as a power-law fit to it: 

Ml(z) = 2.7 x 10 I0 (l+z) 21 . (9) 

This can be compared directly to the res ult obtained from 
the SDSS sample in the local universe by Kauffman n et al.l 
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d2003bl) (compare also iBrinchmann et al.1 [20041) . They find 
that above a stellar mass of 3 x 10 10 M©, the fraction of galax- 
ies with old stellar populations and low star formation rapidly 
increases. Our estimate of the zero redshift brea k mass in the 
star fo rmation rate is M'(0) = 2.7 x 10 10 M Q . iBundy et aD 
(2006) find evidence for an evolving mass limit for star form- 
ing galaxies in the DEEP2 redshift survey at 0.4 < z < 1.4. 
Their best-fitting parametrization of the "quenching mass", 
Mq(z) oc (1+z) 3 ' 5 , is probably comparable to ours given the 
difference in methodology and redshift range, although the 
exponent we find is somewhat smaller. 

In the lower panels of Fig. [3] we plot the star formation 
rate against redshift for different fixed stellar masses. From 
this plot it is apparent that while massive galaxies have the 
highest (absolute) star formation rates, their star formation 
rates drop earlier and quicker than for lower stellar mass ob- 
jects. The star formation rate in logM* ^11 objects drops by 
a factor ~ 200 between z ~ 4.5 and z ~ 1, while objects of 
logM* ~ 9.5 decrease their average star formation rate only 
by a factor ~ 5. Most of the latter decrease happens at z < 1. 
At higher redshifts, low mass galaxies have almost constant 
(average) star formation rates. The more rapid decline we 
find in the star formation rate of massive gala xies is in con- 
trast with the analysis by IZheng et~aT1 d2007l) who combine 
24 fim data with UV data in the COMBO- 17 survey and con- 
clude that the star formation rate declines independently of 
galaxy mass at z < 1 (using stacked detections to account for 
the average star formation rate in sources individually unde- 
tected at 24 /jm). This difference could be due to the (SED- 
based dust corrected) UV light underestimating the true star 
formation rate in logM, > 1 1 galaxies at z < 1 more so than 
at lower masses. Another possibility is some uncertainty in 
the infrared SEDs of these galaxies and the calibration of the 
observed-frame 24 /im-based star formation rate at interme- 
diate redshifts. 

Another fact to note is that the exponent of the star for- 
mation rate as a function of stellar mass at masses below the 
break mass is remarkably stable. In all redshift bins, the ex- 
ponent is close to j3 ~ 0.6, possibly steepening slightly to 
j3 ~ 0.7 at z < 0.5. This also means - since the exponent 
is always smaller than one - that the specific star formation 
(star formation rate per unit stellar mass) rate is a monoton- 
ically decreasing function of stellar mass at all masses and 
redshifts. 

How does this compare t o other studies of star for mation 
rates? In the local universe, IBrinchmann et al.l (|2004) find a 
very similar slope of the star formation rate with stellar mass 
as well (logM* oc 0.6 logM*). 

At intermediate redshift, in fact, iNoeske et al.l (l2007bl) 
argue for the existence of a "main sequence" of star 
forming galaxies that obeys the relation logM* = (0.67 ± 
0.08) logM*- (6.19 ±0.78) at 0.2 < z < 0.7 and 10 < 
logM*/M Q < 1 1 using DEEP2 survey data combining 24 fim 
and optical emission line fluxes. This is in very good agree- 
ment with our results. They also see indications of slight flat- 
tening towards higher redshift, although this cannot be said 
w ith confidence f r om th eir sample. 

Pap ovich et al.l (120061) studied massive galaxies at 1 < z < 
3.5 detected by the Spitzer Space Telescope at 3-24 /im, by 
HST, and in ground based observations in the GOODS field. 
They also find a relation between stellar mass and star forma- 
tion rate, and although they do not attempt to quantify this re- 
lation, their data is consistent with the relation proposed here. 
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FIG. 5. — Star formation histories of galaxies with mass between 10 8 and 
10" 5 Mq, calculated by integrating the star formation rate as a function of 
mass and time (Fig. [5). . 

They also compare their high-z results to similarly measured 
star formation rates in z < 1 galaxies selected from COMBO- 
17 and find similar results to the evolution presented here. 

Also, our star formation rates (and stellar masses) in the Ly- 
man break galaxy regime agre e very well with detaile d studies 
of this class of objects. As in[Papovich et al.l (1200 ll) . we find 
stellar masses of 9.5 < logM* < 10.5 and star formation rates 
(at z ~ 3) between 10 and 100 M Q yr _1 . 

We are therefore confident that our relationship between 
stellar mass and star formation rate is robust and likely to be 
valid in spite of its rather simplistic footing. It is understood 
that there is an underlying broad distribution of star forma- 
tion rates and much intrinsic scatter around this relation, es- 
pecially at higher redshift where one is likely to encounter 
massive starbursts. Also, the relation as shown here is only 
strictly valid for "star forming galaxies", as quiescent objects 
(most notably at lower masses) are most likely to be unde- 
tected (however, they will be missing from both 0(M«) and 
M*{M*) in the present study). 

Finally, given the average star formation rate field 
M*(M*,t), we can ask what the star formation histories look 
like on average as a function of mass, i.e. we can ask for tracks 
along M* and t that solve the equations 

dM»/df = M^M^t) 
M*(r ) = M . 

Figure [5] shows these solutions for masses Mq at to of 8 < 
logM* < 11.5. Note that Mq and fo are defined as the mass 
and time at z = 4.5. In interpreting these as star formation his- 
tories of galaxies, we assume that the dominant contribution 
to the growth of stellar mass is from star formation, i.e. we 
neglect growth by accretion and merging, which we will see 
is a reasonable approximation in the following sections. 

It has long been thought that star formation histories of 
galaxies are roughly exponentials on average, with later type 
and lower mass galaxies having longer e-folding times (e.g. 
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FIG. 6. — The partial derivative with respect to time of the stellar mass function as measured from the parametrization of its redshift evolution (see text). The 
left hand panel shows the derivative as a function of stellar mass at different redshifts; thicker lines indicate higher z with values between z = 0.5 and z = 4.5 
as in Fig.Q] The right hand panel shows the same functions but plotted against redshift, with higher stellar mass indicated by thicker lines; the mass values are 
indicated and the same as in Fig. [3] At masses below M* , where the mass function is a power law, the evolution becomes independent of mass. Therefore curves 
for masses below logM* ~ 10 in the right panel lie on top of each other. 



Sandage 1986). From investigations of the present day stel- 
lar population in ell iptical galaxies, this was confirmed by 
Thom as et alJ ([2005), at least for spheroidals. iHeavens et al.1 
(2004) use spectra of galaxies from the SDSS to infer star for- 
mation histories of galaxies as a function of mass and find 
similar trends. Here, Fig. [5]- from measurements of star for- 
mation rates as a function of stellar mass and redshift - di- 
rectly shows that a relationship between star formation his- 
tory and galaxy mass holds. Rise time, peak star forma- 
tion rate, peak time, and post-maximum (exponential) decay 
timescale are all correlated with mass. More massive galaxies 
have steeper and earlier onsets of star formation, higher peak 
star formation rate, and the following exponential decay has 
a shorter e-folding time. Less massive galaxies have a more 
gradual onset of star formation, a broader peak of star forma- 
tion activity, and their maximum star formation rate is lower 
and occurs later. Their subsequent exponential decline has a 
longer e-folding time. 

Since the star formation rate as a function of mass is well 
represented by a power law in mass with constant slope, we 
can write the early star formation rate (at masses below the 
break mass) as AM /At =A(f)M 13 , with A(t) ocM° fromeq. (O, 
and monotonically decreasing with time. Hence, galaxies 
which begin forming stars later in time, have longer rise times, 
reach lower peak star formation rate, and also lower final 
mass. The latter statement can be derived from the fact that 
the break mass decreases with time, and so galaxies that start 
forming stars later have their star formation quenched at lower 
mass than galaxies that start forming stars early, when the 
break mass is still high. 

As a side note, the existence of the curious universal power- 
law relation between star formation rate and stellar mass 
can be explained using disk galaxy sca ling relations and the 
Schmidt-Kennicutt star formation law dKenn icutt 19983). If 
one takes the Schmidt-Kennicutt law to read Ssfr oc S /id yn 
(with £ denoting surface density and t being the dynami- 



cal timescale, see ISilklll997t lKennicutt)[l998bl) . then M* oc 
R 2 Ssfr oc R 2 EV/R oc E Vfl (with V denoting the character- 
istic velocity of the system and R its radius). Using the virial 
relation we have M» oc M®- 5 R°- 5 , which, with M oc R 3 , be- 
comes M, oc M° 62 . Using the disk scaling relations given in 
iMo et al.1 d!998l) instead of the (halo) virial relation (M oc V c 3 ; 
Rd oc Vc), one obtains M* oc M* for average halo spin param- 
eter and fraction of disk angular momentum. Both of these 
arguments assume a fixed scaling of cold gas surface density 
with total mass surface density. 

5. THE TIME DERIVATIVE OF THE MASS FUNCTION 

Now that we have assembled all the necessary data, we turn 
back to eq. (f2]). The right hand side of eq. (f2]) has two terms. 
The first is the observed time derivative of the mass function, 
the second one is the change with time induced in the mass 
function due to star formation. Here we will deal with the 
first term. 

We use the parametrization of the evolution of the stellar 
mass function as established in § [3] From this parametriza- 
tion, it is a straight forward procedure to calculate the time 
derivative d<fi/dt. For convenience, we will henceforth al- 
ways measure time derivatives in terms of relative changes, 
(1 /4>)(d(f>/dt). These are plotted against mass in the left hand 
panel of Fig.[6l we use thicker lines to indicate higher redshift. 
The right hand panel shows the same quantity, this time plot- 
ted against redshift, with line thickness indicating increasing 
stellar mass. We plot data for the same redshifts and masses 
as in Fig. [3] 

The shaded regions indicate the \-a confidence limits cal- 
culated using Monte Carlo simulations of the fit of the evo- 
lution of the mass function, eqs. (5) - (7) and Fig. [2] in §[3] 
It is important to note that while the uncertainties on the 
mass function itself are rather large, the parameterization of 
its evolution provides much smaller formal uncertainties on 
the derivative with respect to time that is of interest here. 
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FIG. 7. — The relative change per Gyr in thegalaxy stellar mass function due to star formation only, plotted against mass (left panel) and redshift (right panel). 
These curves are the result of evaluating eq.[T] Confidence limits are derived by directly propagating the confidence limits on the average star formation rate 
through eq.ff] 



Note that our parametrization of the mass function assumes 
a fixed faint-end slope of a = —1.3 (the data do not allow us 
to measure the slope at higher redshift; up to z ~ 2, the slope 
is consistent with being constant). Hence the time derivative 
of the mass function becomes independent of mass at masses 
below the characteristic mass, M*(z), where the mass function 
is a power law. In the right hand panel of Fig. [6] this means 
that the low mass curves (below logM* ~ 10) all lie on top of 
each other. 

The most rapid buildup of the mass function happens (not 
unexpectedly) at high redshift and at the high-mass end. At 
Z ~ 4 and logM, ~ 11, the number density approximately 
doubles every Gyr, while the density of logM, ~ 10 objects 
(typical Lyman break galaxies) increases by 50% per Gyr. At 
z < 1, the change in number densities is much more modest, 
of the order of 10-20% per Gyr (and likely decreasing further 
rapidly), so that the mass density roughly doubles between 
Z ~ 1 and z = 0, a behavior observed and confirmed by many 
studies mentioned in the introduction. 

There is rather little differential evolution in the Fors Deep 
Field stellar mass function in the sense that the ratio of the 
change at high mass and at low mass evolves from a factor 
of ~ 4 a t high z to a facto r ~ 2 at z ~ 0.5. This has been 
noted in iDrory et all (120051) : see their Fig. 4. iFontana et al.1 
(2006), analyzing data similarly in the GOODS-MUSIC sur- 
vey, find that this ratio tends to zero at late times, such that 
the buildup of the logM* > 1 1 galaxy population is complete 
by z ~ 1 while it still continues slowly in our sample (com- 
pare a lso studies of elliptical galaxies, e.g. Ivan der Wei et al.l 
120051) . This discrepancy, apparent in a number of datasets - 
whether it is due to the small volume of the FDF at low z, the 
photometric redshifts, the stellar population modeling tech- 
nique, or other effects - has yet to be resolved. 

Since the main focus of this paper is to discuss the method 
presented here to analyze mass function and star formation 
rate data, we will not spend too much time discussing pecu- 
liarities of the present dataset. We are fully aware that the 
current data have significant deficiencies (both in mass and 



star formation rate), yet they still offer the possibility to dis- 
cover and learn about fundamental trends in the formation and 
evolution of galaxies. 

6. THE EVOLUTION OF THE MASS FUNCTION DUE 
TO STAR FORMATION 

After the derivative of the mass function, we turn to the 
second term on the right hand side of eq. [2] We proceed to 
analyze the relative change in the stellar mass function that 



which is 
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is due to star formation only, "0 = 4; ' f) 

times eq. Q] We will use the symbol ip to refer to this quantity 
in what follows. 

To compute the function ip = <fi~ l (— d/dM*) (</>M*), in each 
time (redshift) bin, we use the parameterized fits to </>(M* , t ) 
andM >t (M*,f) listed in Tables Q] and |2j respectively. We plot 
the resulting curves in Fig. [7] The left hand panel shows the 
relative change of the stellar mass function plotted against 
stellar mass in seven redshift bins spanning 0.5 < z < 4.5; 
as in the previous plots, thicker lines indicate higher red- 
shift. The right hand panel shows the same data but now plot- 
ted against redshift with stellar mass coded as line thickness. 
Masses and redshifts are the same as in Fig. [6] Confidence 
regions are derived directly from the confidence limits on the 
average star formation rates propagated through eq.[T] 

rp(M* , t) is positive at all times. This means that the number 
density of galaxies is predicted to increase with time due to 
the star formation activity of the galaxies. Note that this is not 
a trivial statement: If the number of objects leaving a mass 
bin [M^M* +dM*) in the time interval df due to their star 
formation rate M^{M^ , t) is larger than the number of objects 
entering the same bin from lower mass due to their own star 
formation (see eq. (Q}), the number density at this mass would 
decrease. It is a matter of comparing the derivatives w.r.t. 
mass of the mass function and the star formation rate. 

In the low-mass range, the stellar mass function is well rep- 
resented by a power law with slope a ~ —1.3 (locally), and 
we have assumed that this slope does not evolve with time 
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(which is backed b y the data up to z ~ 2; see § [3] and also 
iFontana eT al. 2006). The star formation rate is also well de- 
scribed by a power law (see § 0J, M* oc M* , with /3 ~ 0.6. 
Eq. (0 times -1 hence reduces to 

1 d 

t/>oc — \M"MP] oc -(a + /3)Mt l - (10) 

This term is positive if a+ (3 < 0. With a ~ -1 .3 and (3 ~ 0.6, 
this condition is satisfied at all times at the faint end. At the 
bright end, the mass function steepens quicker than the star 
formation rate does, and the condition remains satisfied. 

We wish to note here that we are very likely to miss low- 
mass low star formation part of the galaxy population, hence 
we very likely underestimate the slope of the stellar mass 
function, especially at higher redshift. 

At high redshift and high mass, we see that star formation 
leads to a strong increase in number density. The number of 
objects with logM, > 1 1 roughly doubles every Gyr at z ~ 
4. This rate of buildup quickly decreases with time, slowing 
down to an increase by 50% per Gyr by z ~ 2.5. At z < 1.5, 
star formation has dropped to contribute a relative increase in 
number density of ~ 0. 1 Gyr" 1 , and to 0.05 Gyr" 1 at z ~ 0.5. 

At low masses, logM* < 10, the number density very 
rapidly increases (due to star formation) at a rate compara- 
ble to that of high-mass objects. Unlike the high mass objects, 
however, the growth due to star formation persists in low mass 
objects for much longer, dropping to a rate of 0.5 Gyr" 1 at 
Z ~ 1 and 0.2 Gyr" 1 at z ~ 0.5. While the number density of 
massive galaxies is not changing significantly anymore due to 
star formation at low z, that of low-mass galaxies does. 

7. THE ASSEMBLY RATE 

Finally, we have all the ingredients necessary to fully eval- 
uate eq. 0, which we do in Fig. [8] showing our results again 
versus mass in the left panel and against redshift in the right 
panel. We will use the less cumbersome term assembly rate 
for "change in number density due to merging" in what fol- 
lows. We remind the reader to be mindful of the difference 
between this quantity and the merger rate (see §0. 

Confidence regions in Fig. [8] are calculated by means of 
Monte Carlo realizations of data within the uncertainties in 
the mass function evolution and the star formation rates, and 
propagating those through eq. 0. It is worth noting that since 
the constraints on the derivatives of <f>* and logM* with re- 
spect to time are much tighter than the constraints on these 
quantities themselves, the error budget in Fig.[8]is dominated 
by the uncertainty in the average star formation rate. This 
leads to formally small errors at masses and redshifts with 
low star formation rates. 

The curves in Fig. [8] are negative at low mass and posi- 
tive at high mass, with the zero-crossing mass evolving from 
higher masses at high redshift to lower masses at lower red- 
shift. This means that once star formation is accounted for, a 
source of high mass objects and a sink of low mass objects is 
necessary to explain the mass function data. In the most gen- 
eral terms, this is consistent with merging activity building up 
higher mass objects from (many) lower mass ones. In this pic- 
ture, the merger rate decreases with time (the absolute value 
of the curves decreases with time). Also, the mass at which 
the merger rate tends to increase the number density instead 
of decreasing it (the mass above which objects are accreted 
onto and below which objects are destroyed) decreases with 
time. However, an important caveat is that a loss of objects 



could also be due to (faint) objects dropping out of the survey 
as soon as their star formation activity ceases. This can partly 
explain the negative values observed at low masses and the 
trend mentioned above. 

In the following we shall look into the merger rates as a 
function of mass and time in more detail. 

7.1. High redshift and high mass 

At high redshift and high masses, the time-evolution in the 
mass function shows a rapid increase in number density of 
high mass objects (see §0. Star formation accounts for most 
of this increase, leaving a difference consistent with zero at 
logM* ~ 1 1 and a small positive difference of 0.2-0.8 Gyr" 1 
at logM* ~ 1 1 .5 and z > 3. However, our parametrization of 
4>*(z) flattens at high redshifts and over-predicts the normal- 
ization of the mass function in the bin at high z (compare the 
high-z data point in the upper panel of Fig. 0. As a result, 
the total change in the mass function is smaller in the param- 
eterized fit to 4>*(z) compared to the raw data by about 30%. 
Therefore, if we use the Schechter functions shown in Fig.Q] 
directly, which show a steeper increase in cf>* at z > 3, instead 
of the smooth parametrization <fi* oc (1 +z)~ 107 , the curves in 
Fig.[8]move up, and the inferred assembly rate becomes larger 
by the same amount. 

Given this uncertainty, we feel that the star formation rate 
we observe is roughly sufficient to explain the increase in 
number density at high mass and high redshift. As an up- 
per limit, at most 0.5-0.8 effective major mergers per Gyr 
are consistent with the data. This number, though, would be 
enough to transform most of these early high mass objects 
into ellipticals contemporaneously with their major star for- 
mation episode. This is consistent with the picture emerg- 
ing from various observations in the local universe and at 
high redshift which show that most massive ellipticals form 
quickly and early, with more massive galaxies being older, 
and their stars haying formed more rapidly (e.g.lBender et al. 
1 1 9961: iTreu et al]l200 It [van de Ven et al JI20031 ; iThomas et al. 
120051: Ivan der Wei et al.ll2005l:lRenzinill2006l) . 

7.2. High redshift and low mass 

As we have seen, the data point towards small galaxies be- 
ing destroyed and large galaxies being built by mergers, with 
the caveat of selection effects being responsible for parts of 
this trend. 

A plausible interpretation instead of invoking a high de- 
struction rate due to mergers, is based on the fact that low 
mass objects with low star formation rates (and therefore red- 
der colors) are very likely missing from the survey. If the 
slope of the mass function were steeper due to the existence 
of a class of low-mass objects missing from the survey, we'd 
need an even higher destruction rate at low masses, due to the 
increase in numbers and decrease in the average star forma- 
tion rate. However, if the star formation rate as a function 
of mass were steeper, M* oc M^ 8 instead of M* oc M^ 6 either 
due to lower star formation at low masses or higher star for- 
mation at high masses, the difference would be close to zero 
below logM* ~ 1 1 at all times, and little change (~ 0.2 Gyr" 1 ) 
in number density due to merging would be necessary. This 
would imply that at low masses, merging is creating objects 
of any particular mass at the same rate as these objects are 
incorporated into larger structures. 

7.3. Low redshift and high mass 
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FIG. 8. — The assembly rate: the relative change in the galaxy stellar mass function that is caused by anything other than star formation. We interpret this 
quantity as being due to accretion and merging of galaxies. The left hand panel is plotted against mass, the right hand panel against redshift. 



Our current dataset is most reliable (both in terms of mass 
completeness and robustness of the star formation rate esti- 
mate) at redshifts z < 2 and at masses logM* > 10.5. Here, 
the photometric data are deep enough to detect objects over a 
large dynamic range in M*/L, and, also, intergalactic absorp- 
tion (the Lyman break) does not affect the optical bands yet, 
such that the object's spectral energy distributions are well- 
sampled from the rest-frame UV into the optical and near-IR. 

At z < 2.5 the assembly rate is increasing towards higher 
stellar masses and towards lower redshift, mostly due to the 
rapid decline of the star formation rate in this regime, but it is 
always very low. Most of the increase in number density of 
massive galaxies is fully explained by the star formation rates 
observed. 

It is very difficult to measure the high mass end of the mass 
function accurately. These objects are rare, and on the steep 
exponential tail of the mass function any small error in dis- 
tance (due to the use of photometric redshifts) or error in 
M^/L will lead to a large error in the number density. Also, 
the effects of cosmic variance are largest in this mass range 
and may heavily bias results especially in relatively small-area 
surveys as the FDF. Some assurance is provided by the fact 
that the z ~ 1 mass functions of a number of surveys agree rea- 
sonably well. If our measurement of the high mass end of the 
mass function at z < 1 .5 is correct, mergers only contribute to 
the growth of the mass function at logM* > 1 1 and only on a 
low level. A growth of only 0.1 - 0.2 Gyr" 1 is consistent with 
our data. This is equivalent to about 1 effective major merger 
per object in the redshift range 1.5 > z > 0. If, on the other 
hand, there is less evoluti on at the high mass end of the mass 
function, e.g. as found bv lFontana et alj d2006l) . then no merg- 
ing at all would be required at high m ass (compare also th e 
mass functions in Fontana et al. 2004 and Drory et alJl2004l) . 
However, there is strong observational evidence that recent 
mer ging does play a role at least in field elliptical galaxies 
(e.g.lSchweizerill982l:IWIiitmore et alJll997HGoudfrooii et alj 
l200lHGenzeletaUl2001h . " 

Other studies of the merger rate do not help to answer this 
question. On one side. Ivan Dok kum (2005) studied signatures 



of interaction and tidal debris in a sample of nearby red galax- 
ies. With some assumptions on the timescale of the merger 
events, they conclude that the present day accretion rate in 
E/S0 galaxies is AM/M = 0.09 ±0.04 Gyr -1 , a value that 
is fully consistent with our measurement, assuming that red 
sequence galaxies comprise the majority of the lo gM„ > 1 1 
local galaxy population. Also, iBell et all d2006bl) find that 
50-70% of galaxies with logM* > 10.7 have underg one a ma- 
jor merge r since z ~ 1 (se e al so [B ell et al. 2006a). On the 
other side. lLin et"aT] (120041) and lLotz etail ((^006) find a much 
lower merger rate. They argue that in the DEEP2 survey, only 
~ 10% of mass ive galaxies have undergone a major merger 
since z ~ 1.2. iBundy et alJ (|2006t) also argue for a low con- 
tribution of mergers at the h igh mass end. An overv iew of a 
number of studies is given in Kartalte pe et al.1 ([2007). 

7.4. Evolution with time 

As we have mentioned, the general trend in the data is for 
the assembly rate to become more positive with time. The 
zero-crossing mass decreases with time. In high mass objects, 
the assembly rate starts out at zero or slightly below, with 
star formation rates dominating the initial increase in number 
density. The assembly rate then increases, such that at low 
z, merging becomes the dominant process that increases the 
number density at masses logM* > 11. The growth rate in 
this mass and redshift range is 0.1 - 0.2 Gyr" 1 . Lower-mass 
galaxies are consumed in these mergers. The mass at which 
the assembly rate transitions from objects being consumed to 
objects being formed decreases with time, roughly tracing the 
break mass in the star formation rate. The zero-crossing of 
the assembly rate occurs at logM* = 10.8 at z = 0.5, logM* = 
11.0 at z= 1.5, logM* = 11.3 atz = 2, and logM* = 11.5 at 
Z = 3.5. While merging is still happening below this transi- 
tion mass, the net number of objects decreases there. Hence, 
at lower redshift, the number density of progressively lower- 
mass galaxies starts to increase due to merging. This is a di- 
rect observation of the top-down buildup of t he red sequence , 
where major merger remnants are found (e.g.lBell et alj|2004t 
iFaber et alJl2005MBundv et"aT]l2006l:lBell et al.ll2007l) . 



12 



Drory & Alvarez 



It is interesting to note that in CDM, e.g. in the 
iPress & SchechterTd 19741) formalism, we expect the opposite 
behavior: the characteristic mass increases with time, and 
hence progressively more massive halos preferably get in- 
corporated into even larger structures, and their number den- 
sity begins to decline. However, this formalism assumes that 
merged halos do not survive as substructure within their par- 
ent halos. Accounting for the survival of substructure, e.g. 
the existence of satellite galaxies or galaxies within clusters, 
might account for this apparent discrepancy. This will be ex- 
plored in a follow-up paper. 

8. SUMMARY, DISCUSSION, AND CONCLUSIONS 

Summarizing our results, we present average star formation 
rates as a function of stellar mass and redshift, M*(M*,f), for 
galaxies in the Fors Deep Field spanning < z < 5. We di- 
rectly measure their star formation histories as a function of 
stellar mass by integrating the former quantity. 

At z > 3 the average star formation rate is monotonically 
rising with stellar mass and is well-fitted by a power law of 
stellar mass with exponent (3 ~ 0.6. The average star forma- 
tion rates in the most massive objects at this redshift are 100 
- 500 M yr _1 , suggesting that these are the progenitors of 
modern day massive ellipticals seen during their major epoch 
of star formation. At redshift z ~ 3, the star formation rates 
starts to drop at the high mass end while still being an increas- 
ing function of mass at lower masses. The power-law expo- 
nent is consistent with staying constant. As redshift decreases 
further, the star formation rates drop at progressively lower 
masses (downsizing), dropping most rapidly with redshift for 
high mass (logM* > 11) galaxies. The mass at which the 
star formation rate starts to deviate from the power law form 
(break mass) progresses smoothly from logM 1 > 13 at z ~ 5 
to log M 1 ~ 10.9 at z ~ 0.5. Fo r comparison, in the local uni- 
verse [kauffman n et alj d2003bl) find a characteristic mass of 
logM* ~ 10.4 at which the galaxy population changes from 
actively star forming at lower stellar mass to typically quies- 
cent at higher masses. The break mass is evolving with red- 
shift according to M\(z) = 2.7 x 10 10 (1 +z) 21 . 

We directly show that a relationship between star formation 
history and galaxy mass holds. Rise time, peak star forma- 
tion rate, peak time, and post-maximum (exponential) decay 
timescale are all correlated with mass. More massive galax- 
ies have steeper early onsets of star formation, their peak star 
formation occurs earlier and is higher, and the following ex- 
ponential decay has a shorter e-folding time. 

We then present a method to infer effective merger rates (or 
mass accretion rates) from observations of the galaxy stellar 
mass function and the star formation rate as a function of mass 
and time. We quantify the evolution of the mass function due 
to the star formation activity of galaxies. We compare this 
evolution to the observed time derivative of the mass func- 
tion and interpret any difference between the two as due to 
changes in galaxy mass due to external causes; accretion and 
merging are dominant among such processes. We introduce 
and analyze a quantity that we call the assembly rate, which 
is defined as the change in number density of galaxies due to 
accretion and merging. 

We observe sufficient star formation to explain the increase 
in number density at high mass and high redshift. As an upper 
limit, at most 0.8 effective major mergers per Gyr are consis- 
tent with the data. This number, though, would be enough 
to transform most of these high mass objects into ellipticals 
contemporaneously to their major star formation episode. 



The data point towards small galaxies being destroyed and 
large galaxies being built by mergers, with the caveat of selec- 
tion effects being responsible for parts of this trend. A plausi- 
ble interpretation, instead of invoking a high destruction rate 
due to mergers, is based on the fact that low mass objects with 
low star formation rates (and therefore redder colors) are very 
likely missing from the survey. The general trend in the data 
is for the assembly rate to become more positive with time. 
The zero-crossing mass decreases with time. In high mass 
objects, the assembly rate starts out at zero or slightly be- 
low, with star formation rates dominating the initial increase 
in number density. The assembly rate then increases, such 
that at low z, merging becomes the dominant process that in- 
creases the number density at masses logM* > 1 1. A growth 
of 0.1 - 0.2 Gyr" 1 is consistent with our data at z < 1.5. This 
is equivalent to about 1 effective major merger per object in 
the redshift range 1 .5 > z > 0. 

Synthesizing the above into a picture of the formation of 
massive ellipticals, a picture can be drawn in which ellipti- 
cal galaxies would have formed their stars quickly and early 
(z > 2.5) and subsequently quenched their star formation, 
with the highest mass galaxies forming first, and on the short- 
est timescale. Their elliptical morphology may be a conse- 
quence of early merging activity given the merger rate that 
is consistent with our data at high redshift (~ 0.8 Gyr" 1 ; see 
above). Later, these galaxies may continue to evolve along the 
red sequence, however during this phase at z < 1 .5, not much 
growth by merging is happening any more (at most a factor of 
two in stellar mass). 

At more moderate masses of 10 < logM* < 11, we find 
that the evolution of the number density due to mergers points 
towards these galaxies being preferably destroyed at early 
times, while at later times the change in their numbers turns 
positive. At lower redshift, the number density of progres- 
sively lower mass galaxies starts to increase due to merging. 
This is a direct confirmation of the top-down buil dup of the 
red sequence suggested by recent observations (e.g.lBell et alJ 
l2004t iFaber et al] 120051 iBundv et alJl2006l iBeU et alJl2061) 
and the simultaneous shutdown of star formation in these ob- 
jects. 

The results presented here do depend on some assumptions 
that we had to make, mainly due to general deficiencies of 
current observations and more specifically of the dataset we 
have at hand. Assuming that the faint-end slope of the MF 
does not evolve is consistent with the observati ons only up 
to z < 2 dDrorvet all 120051; iFontana et aT]|2006l) . A steeper 
faint-end slope is possible due to faint objects (with lower star 
formation rates than the detected ones) being missed due to 
the survey flux limit. Therefore, a steeper faint-end slope of 
the MF also means a flatter average star formation rate, and 
hence a larger assembly rate at low masses. 

The fitting function one chooses to use to describe the evo- 
lution of the mass function also affects the results. For ex- 
ample, a power-law fit to the evolution of cj>* constrains the 
derivative of the MF w.r.t. time, and thus limits the possible 
range of assembly rates. Ideally, one would like to have MF 
data robust enough to be used to measured derivatives directly 
and avoid the use of fitting formulae entirely. 

The choice of a single average value for the SFR at each 
mass and redshift also affects the results. In reality, galaxies 
show a range of activity at fixed mass. Allowing for such a 
range of possible values instead of a single number would 
lend the method more flexibility, but also demands much 
higher quality star formation data with better-understood se- 
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lection effects. 

The current generation of ongoing large-area surveys of- 
fer the possibility to further explore the build up of the red 
sequence. With a much larger number of red and massive 
objects, one could improve on the analysis presented above 
in various ways. Studying merger processes as a function of 
the environment is one such possibility. Another promising 
approach is to separate galaxies into passive objects and ob- 
jects that are actively forming stars. One can then apply our 
method to these two classes of objects individually, adding 
another term to eq. (O that describes the rate of migration be- 



tween the two sequences. This will allow to investigate the 
relation between merging activity and quenching of star for- 
mation in a very direct way. 



We wish to thank Peter Schuecker, Roberto Saglia, and 
Ralf Bender for stimulating and encouraging discussions. We 
thank Ulrich Hopp for commenting on an early version of the 
manuscript. We thank the anonymous referee for providing 
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